Singular solutions of a fully nonlinear 2x2 system 

of conservation laws 



Henrik Kalisch § Darko Mitrovic §1^ 



Abstract: Existence and admissibility of 5-shock solutions is discussed for the 
non-convex strictly hyperbolic system of equations 

d t u + d x (^)=0 
d t v + d x (v(u - 1)) = 0. 

The system is fully nonlinear, i.e. it is nonlinear with respect to both unknowns, and 
it does not admit the classical Lax-admissible solution for certain Riemann problems. 
By introducing complex-valued corrections in the framework of the weak asymptotic 
method, we show that a compressive 5-shock solution resolves such Riemann prob- 
lems. By letting the approximation parameter tend to zero, the corrections become 
real-valued and the solutions can be seen to fit into the framework of weak singular 
solutions defined in |T2] . Indeed, in this context, we can show that every 2x2 system 
of conservation laws admits 5-shock solutions. 
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1. Introduction 

The main subject of this paper is a system of conservation laws appearing in the 
study of plasmas. The system is known as the Brio system and has the form 



The system is strictly hyperbolic; it is genuinely nonlinear at {(u, v ) : u£ IR, v > 0} 
and {(u, v ) : « 6 1, K 0}, but not on the whole of IR 2 . The system was introduced 
in [2] and thoroughly considered in |15j . There, it was found that for certain initial 
data no solution consisting of the Lax-admissible elementary waves (shock and rar- 
efaction waves) exists. In [15] , Riemann problems for fll.ip were compared to Riemann 
problems for the system 



Numerical computations of appropriate viscous profiles for fll.ip and fll.2j) demon- 
strated surprising similarities. In the same paper, it was shown that certain Riemann 
problems for ( 11. 21) admit 5-shock wave solutions. However, the same fact could not 
be established for any Riemann problem corresponding to (II. ip . In the present work, 
we aim to resolve the question of existence of 5-shock wave solutions of f ll.ip . and the 
question of physical justifiability of such solutions to the Riemann problem associated 
to ( II. ip . We remark that for (II .2p . if the 5 distribution is a part if the solution then 
it is adjoined to the function v (with respect to which the system is linear). However, 
in the case of system (II. ip . our investigation shows that it is more natural for the 5 
distribution to be a part of the function u. 

The study of singular solutions of systems of conservation laws was initiated by 
Korchinski [23] arid Keyfitz and Kranzer [21, 22J. In the last few years, interest in 
the topic has grown, and a sample of results may be found in [31 EJ H31 HSJ HEJ [191 
1201 EMI 1261 1271 1301 1351 1361 138] . One convenient tool for constructing singular solutions 
is the method weak asymptotics. This method has been used recently to understand 
the evolution of nonlinear waves in scalar conservation laws as well as interaction and 
formation of 5-shock waves in the case of a triangular system of conservation laws 
[HI [TOl [12] . We refer the reader to [32] and the references contained therein for further 
applications of the weak asymptotics method. 

In the present work, we introduce an extension of the weak asymptotics method 
to the case where complex-valued corrections are considered for the approximate 
solutions. Even though the imaginary parts of the solutions so constructed vanish 
in an appropriate limit, it appears that considering complex- valued weak asymptotic 
solutions significantly extends the range of possible singular solutions. 



d t u + d x (^)=0, 
d t v + d x (v(u - 1)) = 0. 



(1.1) 



dtu + d x (£)=0 
d t v + d x (v(u - 1)) = 0. 



(1.2) 
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It appears that the weak asymptotic method has so far only been used to construct 
singular solutions of systems for which the flux functions were linear with respect 
to the unknown function which contains 5-distribution. In contrast, note that the 
flux (f(u,v), g(u,v)) = {{u 2 + v 2 )/2, v{u — 1)) associated with the system fll.ip is 
nonlinear in both u and v, and none of the existing methods yield singular solutions 
of this system. Thus it appears that the use of complex-valued corrections is essential 
in the construction of singular solutions for fll.ip . 

Let us next define what we mean by complex-valued weak asymptotic solution, and 
highlight some methods to restrict the notion of solution with the goal of obtaining 
uniqueness. First we define a vanishing family of distributions. 

Definition 1.1. Let f £ (x) G be a family of distributions depending on e G 

(0, 1), We say that f £ = op/(l) if for any test function <f>(x) G T>(M.), the estimate 



The estimate on the right-hand side is understood in the usual Landau sense. Thus 
we may say that a family of distributions approach zero in the sense defined above if 
for a given test function <ft, the pairing (f £ , <ft) converges to zero as e approaches zero. 
For families of distributions f £ (x,t), we write f £ = oxv(l) C T>'(M) if the estimate 
above holds uniformly in t. More succinctly, we require that 



where the function g depends on the test function <p(x, t) and tends to zero as e — > 0, 
and where Ct is a constant depending only on T. We define weak asymptotic solutions 
to a general system of two conservation laws 



Definition 1.2. We say that the families of smooth complex-valued distributions {u £ ) 
and (v £ ) represent a weak asymptotic solution to fll.3p if there exist real-valued dis- 
tributions u, v G C(R + ; V{R)), such that for every fixed t G M + 



(f £ ,(f>) = o(l), as £^0 



holds. 



(f £ (-,t)^)<C T g(e)foite[0,T] 



d t u + d x f(u,v) =0, 
d t v + d x g(u, v) =0, 



(1.3) 



as follows. 




v 



as £ — > 



in the sense of distributions in V'{R), and 



d t u £ + d x f(u £ ,v £ ) 
d t v £ + d x g(u £ ,v £ ) 



ot»(1), 

Op/ (1). 



(1.4) 



4 



It is evident that this definition requires some additional assumptions of the fluxes 
/ and g. In particular, / and g must have an extension into the complex plane. 
One may for instance restrict to fluxes that are real- analytic, though in principle a 
wider class of fluxes is possible. The main issue in the requirement on the fluxes, and 
indeed with this method of constructing solutions is the question of uniqueness. For 
example, by adding a constant term of order 0(e) to any weak asymptotic solution, 
one immediately obtains two different weak asymptotic solutions which correspond 
to the same solution if a more restrictive concept is used. 

One way to narrow the class of solution candidates is to require distributional 
solutions to satisfy the equations in a stronger sense than the one defined in Definition 
11.21 This approach entails substituting them into ( jl.ip . and to check directly whether 
the equations are satisfied. This strategy involves the problem of multiplication of 
singular distributions. The problem of taking products of singular distributions was 
overcome by Danilov and Shelkovich in [TJ] in a rather elegant way. In their work, 
the weak asymptotic solution is constructed such that the terms that do not have a 
distributional limit cancel in the limit as e approaches zero. As a result, it is not 
necessary to include singular terms in the definition of the weak solution. Thus, the 
problem of multiplication of distributions is automatically eliminated, and the class 
of possible solution is significantly reduced. 

There are also several other reasonable ways to multiply Heaviside and Dirac dis- 
tributions. In [T7J [37], a number of definitions of weak solutions of (jl.3p are 
introduced. Among the latter approaches, we emphasize the measure-type solution 
concept introduced in [7_J H7J . Moreover, the framework from [T7j yields uniqueness of 
solutions if an additional condition of Oleinik-type is required, and that is probably 
the only work so far which obtains a uniqueness result for arbitrary initial data in a 
class of distributional solutions weak enough to allow delta distributions. However, 
uniqueness has also been obtained for special classes of initial data by LeFloch in [25] , 
and by Nedeljkov [30] . 

We remark that a systematic study of multiplication of distributions problem is 
investigated in the Colombeau algebra framework [51 [HI [28] . In these works, problems 
of the type considered here are also investigated. Actually, Definition 11.21 can be 
understood as a variant of appropriate definitions in [6j [29J [31] . The main difference 
is that in the present solution is found pointwise with respect to t e M+, 

and it is required that the distributional limit of the weak asymptotic solution be a 
distribution. The latter is not necessary in the framework of the Colombeau algebra 
though it may be tacitly assumed. 

The plan of the present paper is as follows. We will provide a review of the definition 
of weak singular solutions from [12] in Section 2. It turns out that a somewhat 
more general statement is appropriate here. Moreover, it will be proved that any 
2x2 system of hyperbolic conservation admits singular solutions of this type. In 
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Section 3, weak asymptotic solutions of the Brio system are found. The results of 
that section are very important since they represent a justification of the concept 
introduced in Section 2 which will be applied in the Section 4. In that section, it is 
shown that the limit of the weak asymptotic solutions satisfy the equation in the sense 
of Definition 12.11 Also, an adaptation of the Lax admissibility concept is proposed 
which provide physically sustainable solutions to corresponding Riemann problems. 
The final section is the Appendix where we consider other possibilities for existence 
of 5-shock solutions. 



2. Generalized weak solutions 

In this section, the definition of weak singular solutions of a 2 x 2 system of conser- 
vation laws provided in [12] is reviewed. Indeed, we shall show that any 2x2 systems 
of the form 

d t u + d x f(u,v) =0, 
d t v + d x g(u,v) =0, 

admits <5-type solution in the framework introduced in |12j . While the definition in 
[12] is given only for solutions singular in the second variable, while assuming that the 
flux functions / and g are linear in the second variable, it appears that the definition 
can actually be made more general. Suppose r = {7^ | i G /} is a graph in the closed 
upper half plane, containing Lipschitz continuous arcs 7$, i G I, where J is a finite 
index set. Let Iq be the subset of I containing all indices of arcs that connect to the 
a;- axis, and let r = {x° k \ k G Jo} be the set of initial points of the arcs 7^ with 
k G Jo- Define the singular part by a(x, t)S(T) = J2iei a i( x i t)^(li)- Let (u,v) be a 
pair of distributions, where v is represented in the form 

v (x, t) = V(x, t) + a(x, t)5(T), 

and where w, V G L°°(M x M+). Finally, the expression denotes the tangential 

derivative of a function ip on the graph 7$, and connotes the line integral over the 
arc 7j. 

Definition 2.1. a) The pair of distributions u and v = V + a(x,t)5(T) are called 
a generalized 5 -shock wave solution of system (11. 3p with the initial data Uq{x) and 
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Vo(x) + Y^j «fc(^o; 0)5 (x — x°) if the integral identities 

(ud t ip + f(u,V)d x <p) dxdt+ / U (x)<p(x, 0) dx = 0, (2.1) 





(^%> + 0(u, I03r¥>) ^xtft (2.2) 

+ ]T /" ai(ar,*)^ + / V (x)<p(x,0) dx + ^ a fe (x°, 0)<p(x° k , 0) = 0, 
ie/ ^ ^ R fce/ 

/ioW /or all test functions ip G X>(R x R+). 

The next definition concerns the similar situation where the singular solution is 
contained in u, and v is a regular distribution. Thus we assume the representation 

u(x, t) = U (x, t) + a(x, t)5(T), 

where now U,v G L°°(IR x K + ), and a(x, t)5(r) is defined as before. 

Definition 2.1. b) The pair of distributions u = U+oi(x,t)6(T) andv is a generalized 
5-shock wave solution of system (\1.3\i with the initial data Uq{x) + afc(xg, 0)5 (x — 
x®) and Vo(x)) if the integral identities 

(Udt<p + f(U,v)d x (p) dxdt (2.3) 
+ J2 [ a l (x,t)^^+ [ U (x)ip(x,0)dx + J2M4,0Mx° k ,0) = 0, 

{yd t ip + g(U,v)d x <p) dxdt + / V (x)ip(x,0) dx = 0, (2.4) 





hold for all test functions <p G X?(IR x M+). 

This definition may be interpreted as an extension of the classical weak solution 
concept. Moreover, as noticed in e.g. [TJ, the definition is consistent with the concept 
of measure solutions [7J [T7] . 

Definition 12. II is quite general, allowing a combination of initial steps and delta dis- 
tributions; but its effectiveness is already demonstrated by considering the Riemann 
problem with a single jump. Indeed, for this configuration it can be shown that a 
5-shock wave solution exists for any 2x2 system of conservation laws. Consider the 
Riemann problem for (II. 3ft with initial data u(x, 0) = Uq(x) and v(x, 0) = Vq(x), 
where 

tt i \ x<0 \r i \ l Vl > x <0 roti\ 

U {x) = < , V (x) = < (2.5) 

I «2, X > \V2, X > 

Then, the following theorem holds. 



Theorem 2.1. a) Ifu\ 7^ u 2 then the pair of distributions 

u(x,t) = U (x-ct), (2.6) 
v(x,t) = V (x-ct) + a(t)5(x-ct), (2.7) 

where 

c = [f{ "'p ] = ^ 2 )-/K^) and a{t) = {c[V ] - [g(u, V)])t, 

[u] U2-U1 

represents the S-shock wave solution of ( 11.31) with initial data Uq(x) and Vq(x) in the 
sense of Definition ^. 1\ a). 

Theorem 2.1. b) If Vi ^ v 2 then the pair of distributions 

u(x,t) = U (x-ct) + a(t)5(x-ct), (2.8) 
v(x,t) = V (x-ct), (2.9) 

where 

[g{U,v)\ g(u 2 ,v 2 )-g(u 1 ,v 1 ) , , , r n r / , n 
c= — n — = , a{t) = {c[U]-[f{U,v)])t 

represents the 8-shock solution of (I1.3P with initial data Uq(x) and Vq(x) in the sense 
of Definition \2.1\ b). 

Proof. We will prove only the first part of the theorem as the second part can be 
proved analogously. We immediately see that u and v given by (12. 6p and (12. 7p satisfy 
( 12. ip since c is given exactly by the Rankine-Hugoniot condition derived from that 
system. By substituting u and v into (12. 2p . we get after standard transformations: 

(-c[V] + [g(u, V)]) (p(ct, t) dt- I a'(t)(p(ct, t) dt = 0. 

Jk + 

From here and since a(0) = 0, the conclusion follows immediately. □ 

As the solution framework of Definition 12.11 is very weak, one might expect non- 
uniqueness issues to arise. This is indeed the case, and the proof of the following 
proposition is an easy exercise. 

Proposition 2.1. System (II. 3ft with the zero initial data: u\ t= o = i>|t=o = admits 
S-shock solutions of the form: 

u(x, t) = 0, v(x, t) = (35(x — C]t) — (35(x — c 2 t), 

for arbitrary constants (3, c\ and c 2 . 

At the moment, we do not have a general concept for resolving such and similar 
non-uniqueness issues. In the case of the Brio system which we shall consider in the 
sequel, we are also not able to obtain uniqueness, but we can prove that there always 
exists a physically reasonable solution to the corresponding Riemann problem. 
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Finally, let us remark, that it is of course possible and reasonable to give a definition 
along the lines of Definition 12.11 which allows for simultaneous concentration effects 
in both unknowns u and v. In this generalized 5-shock wave solution of (jl.3p 

with the initial data U (x) + J2i a k{ x oi 0)5[x — x®) and V (x) + J2i 0k( x o-> 0)<5(x — x®) 
would have the form u = U + a(x, t)5(T) and v — V + (3(x, t)5(T), and satisfy 

(Ud t ip + f(U,V)d x if) dxdt (2.10) 

+ Y)/ a l (x,t)^^ + [ U (x)<p(x,0) dx + J2a k (x o k ,0)<p(x o k ,0) = 0, 
iei J ^ Jr kei 

(Vd t ip + g(U,V)d x( p) dxdt (2.11) 





+ J2 [ ft(M)^ 1 + / V (x)<p(x,Q) dx+J2^(xl0Mxl0) = 0, 

te 

in [I]. 



for all test functions ip G X>(R x M. + ). An example of such a situation can be found 



3. Weak asymptotics for the Brio system 

In this section, we shall construct weak asymptotic solutions for the Riemann prob- 
lem associated to the Brio system (jl.ip . and then show that the weak asymptotic 
solution converges to the generalized weak solution to the system in the sense of Def- 
inition 12.11 This construction is very important since the fact that it is possible to 
find a sequence of smooth approximating solutions to (11. ip . (12. 5p converging to the 
5-shock solution represents a justification of the concept laid down in Section 2. In 
particular, observe that the vanishing viscosity approximation is a special case of the 
weak asymptotic approximation since the term clearly of order Ow(e). 

To find the weak asymptotic solutions we need to find families of smooth functions 
(u £ ), (v e ), such that 

d t v £ + d x {v £ {u £ - 1)) = 025/(1), 

u £ — ^ u, v £ — v as e — > 0, (3.2) 

and such that u(x, 0) = Uq(x) and v(x, 0) = Vq(x) are given by (12. 5p . We shall prove 
the following theorem. 

Theorem 3.1. a) If u\ ^ U2 then there exist weak asymptotic solutions (u E ), (v £ ) of 
the Brio system (ll.ip . such that the families (u £ ) and (v E ) have distributional limits 

u(x,t) = U (x — ct), (3.3) 
v(x,t) = V (x - ct) + a(t)5(x - ct), (3.4) 



where 

c= — r and a(t) = - (c(w 2 - fi) + - 1) - f 2 (W2 - I)))*- 

2(«i - m 2 ) 2 

(3.5) 

Theorem 3.1. b) If vi ^ v 2 then there exist weak asymptotic solutions (u e ), (v £ ) of 
the Brio system (II. ip . sitc/i that the families (u £ ) and (v £ ) have distributional limits 

u(x,t) = U (x — ct) + a(t)6(x — ct), (3.6) 

v{x,t) = V (x-ct), (3.7) 

where 

c= Mui-l)-V2(u 2 -l) ^ a{t)= ( c{u2 _ Ul) + *l + *l-*Z-* 2 2 ) t (3 .8) 

vi-v 2 V 2 / 

Proof, a) Let p G C^°(M) be an even, non-negative, smooth, compactly supported 
function such that 



suppp C (— 1, 1), / p(z)dz — 1, p > 0. 

We take: 

i? £ (x,t) = ip((s-ct-2e)/e)-ip((a;-ct + 2e)/e), 
£ e (x,t) = ip((x - ct - 4e)/e) + \p({x -ct + 4e)/e). 
Next, define smooth functions and VC such that 



(3.9) 



Wi, x < ct — 20s, 
U £ (x, t) = { c + 1, ct - lOe < x < ct + lOe, 
w 2 , x > ct + 2(k, 

I v\, x < ct — 20e, 
K(x, t) = < 0, ct - lOe < x < ct + lOe, 
lf 2 , x>ct + 20s. 

Notice that 

R £ -± 0, U £ R £ -± 0, and U £ 5 £ -± 2(c + l)S(x - ct). (3.10) 
Moreover, we have 

V £ R £ = 0, V e 8 £ = 0, and 5 £ R £ = 0. (3.11) 

Now make the ansatz 

u £ (x,t) = U E (x,t), 

v £ (x,t) = V £ (x,t) + a(t)(6 £ (x,t) + R £ (x,t)), 



(3.12) 
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and substitute it into equations (I3.ip . Notice first of all that 

v 2 (x,t) = V 2 + a 2 (t)(R 2 £ +5 2 ) 

by invoking ( 13. lip . Focusing on the expression R 2 £ + S 2 , we take ip e C C °°(M) and 
consider the integral 

(R 2 £ + 5 2 )<p dx 

i ( - P 2 ((x - ct + 2e)/e) - p 2 ((x - ct - 2e)/e) 

+ p 2 ((x - ct + Ae)/e) + p((x - ct - 4e)/e)j<p dx = 0{e). 
In the above reasoning, use was made of the following computation. 

■p (p 2 ((x - ct + ae)/e) + p 2 ((x - ct - fie)/e)) ip(x) dz 
\p 2 {z) (cp(ct + e{z - a)) + <p(ct + e{z + /?))) dz 
±p 2 (z) (2<p(ct) + eip'{ct){p - a)) dz + 0(e), for a,/3 ER. 



The last relation was found by making the changes of variables (x — ct + ae)/e = z 
and (x — ct — f3e)/e = z, and observing that j zp 2 (z)dz = since p is an even function. 
In the case at hand, we use a = (3 = 2 for the first integral, and a = (3 = 4 in the 
second integral. Finally, it becomes plain that 

v 2 = V 2 + o v ,(l). (3.13) 

Therefore, taking into account Definition 11.11 from the first equation in (13. ip . we 
conclude that we need to check whether 

U 2 + V 2 
d t U £ + d x £ Z e = ov.il), 



and this reduces to 

/ (-c[U} + l[U 2 + V 2 ])p(ct,t)dt = o(l), (3.14) 
Jo 

where [U] = u 2 — U\ and [U 2 + V 2 ] = u\ + v\ — u\ — u\. 

However, this is indeed satisfied thanks to the choice of the constant c which was 
found from the Rankine-Hugoniot condition for the first equation in (II. ip . 

Let us now consider the second equation in (13.11) . First, notice that 

d x (v e (u e - 1)) = 8 X (U e V £ + (c + l)a(t)5 E -V e - a(t)5 £ ) + 0^,(1) 

= (fi(l — Ui) + v 2 (u2 — l)6(x — ct) + ca(t)S'(x — ct)) + op'(l). 



11 

Next, note also that 

d t v £ = —c{y 2 — Vi)8(x — ct) + a'(t)5(x — ct) — ca(t)S'(x — ct) + oth(1). 
Adding the latter two expressions, we obtain 
d t v £ + d x {v £ (u £ - 1)) 

= (— c(v 2 - vi) + a'(t) + (vi(l - Mi) + v 2 (u 2 - 1)) S(x - ct) + oxv(l). 

From here, we conclude that choosing a as given in (I3.5p . the first equation in (13. ip 
is satisfied, as well. This concludes the proof of part (a). 

b) In this case, an appropriate weak asymptotic solution is given by 

u £ (x,t) = U £ (x,t) + a(t)(5 £ (x,t) + R l£ (x,t)) + y/2ca{t)Ru(x,t)), 
v £ (x,t) = V £ (x,t), 

where 

VxU! - V 2 U 2 ( ul + vf-uj-l 

— 1 and a(t) = I c[Ui — u 2 ) 



Vi - v 2 



and 



i, x < ct — 20e 
U £ (x, t) = { 0, ct - lOe < x < ct + lOe 
«2, x > ct + 20s 

1, x < ct — 20e 
V £ (x, t) = {0, ct - lOe < x < ct + lOe ; 

2 , x > ct + 20e 

Ri £ (x, t) = -p((x — ct — 2s) /e) — -p((x — ct + 2e)/e); 
s s 

R 2£ (x,t) = — [p((x - ct)/s)}^ ; 

V s 

5 E (x,t) = -p{(x - ct- As)/s) + -p((x -ct + 4e)/e), 

where p is the same smooth non-negative even function as used in the previous ex- 
amples. The proof then follows the ideas of the proof of (a). □ 

An important corollary (to be used in the Appendix) of the proof of the previous 
theorem is that it gives another interesting class of weak asymptotic solutions to (II. ip 
having the 5 distribution as their limit. 

Corollary 1. If U\ = u 2 and v\ = v 2 then for any ceK the families (u £ ) and (v £ ) 
given by (I3.12p are the weak asymptotic solution to (ll.ip . (12. 5p . 
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Proof. It is enough to notice that (13.141) is satisfied independently on c since [U] = 



To close the section, we should mention that while the extension of the weak asymp- 
totics method to complex-valued solutions was crucial for finding a solution of the 
system ( 11. ip . it might not be appropriate in other contexts as it might lead to strong 
non-uniqueness. For example, using complex-valued weak asymptotic solutions of 
similar form as ( 13. 6 P for the inviscid Burgers equation, one may construct a family of 
distinct solutions emanating from the same initial data, all of which also satisfy the 
Lax admissibility condition. 

4. Generalized weak solutions for the Brio system and the 

uniqueness issue 

By comparing Theorem 12.11 and Theorem 13.11 we see that the limit distributions 
u and v given in Theorem 13.11 represent 5-shock solutions to ( II. ip with initial data 
u(x, 0) = Uq(x) and v(x, 0) = Vo(x). However, we want to incorporate such solutions 
into the Lax admissibility concept and this is the goal in this section. 

We focus on Definition 12.11 b) from [12] where the fluxes / and g are given by 
the Brio system. The same can be done with Definition 12.11 a) but it appears that 
the solutions which it generates do not fit into the Lax admissibility concept. More 
details of this case will be provided in the Appendix. 

Recall that in the case v% > > v 2 there exists no Lax-admissible solution to the 
Riemann problem ( II. ip . with the Riemann initial data Uq and Vq given by (12. 5 p (see 
[15]). If V\ and v% do not satisfy this relation, we have the classical Lax admissible 
solution to the appropriate Riemann problem consisting of the elementary waves, i.e. 
shock and rarefaction waves. For L°°-small data, such solution is unique since the 
system is genuinely nonlinear for v > and v < 0. Theorem 12.11 states that we can 
also have 5-shock wave solutions, but as Proposition 12.11 shows, there is strong non- 
uniqueness. In order to eliminate at least some of solutions which are inconsistent 
with the physical intuition, we shall use the Lax compressivity conditions for the 8- 
shock wave. In order to introduce them, let us recall that the characteristic velocities 
for the Brio system are [15] : 



[U 2 + V 2 



] = o. 



□ 



Ai(«, v) = u - 1/2 - ^/1/4 + v 2 , X 2 (u, v) = u - 1/2 + y/l/A + v 2 . 



The corresponding rarefaction waves are given by 



RW 2 : u = - 



RWt : u = -- 




2: 
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The shock waves are given by 

SW h2 : u-u 1 = (l =F y/iv + vtf + l) . 

v + v i V / 

A phase space picture for a given left and right state is shown in Figure 1. 




Figure 1. L is the left, and R is the right state. 



The following definition introduces a compressivity demand on the characteris- 
tics of (II. ip meaning that the characteristics enters the 5-shock from both sides. It 
is standard for the classical shock waves and they are known as Lax admissibility 
conditions. Note the usual demand on the 5-shock wave is an overcompressivity 
condition demanding that both characteristic field Ai and A2, satisfy (14. ip below 
[TO] H~2l [TBI I2IT [36] . However, we were not able to find solutions involving overcom- 
pressive 5-shocks and we confine ourselves on a less restrictive demand which still 
includes concentration effects. The definition concerning the admissible 5-shock so- 
lutions of (II. ip such as defined in Theorem 3.1 follows. 

Definition 4.1. A 5-shock solution of (II. ip . connecting a left state L = (ui,vi) and 
a right state R = (u2,V2) is i- admissible if 

Ai(u 2 , v 2 ) < c < Xi(u!, vx), (4.1) 

for i = 1 or i = 2. For such S shock wave we say that it is compressive. 

Thus for a general Riemann problem, one may say that a solution of (II. ip . (12. 5p 
which contains a 5-shock wave is admissible if it consists of a combination of the 
classical Lax admissible simple waves (shock or rarefaction) and compressive 5 waves. 
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The following lemma will be crucial for existence of admissible 5-shock solutions to 
Riemann problems corresponding to ( 11. II) . 

Lemma 4.1. Assume that the initial data (12.51) are such that U\ = u 2 = u, v\ = 

and v 2 < 0. Then, the S -shock solution 

u(x, t) =u + a(t)5(x — ct), 

(4-2) 

v(x, t) =0, 

where a(t) and c are given by ( 13.81) . represents 1-admissible S -shock solution. 

Proof. The functions given by (14. 2 j) represent <5 shock solution to (11.11) . (12.51) according 
to Theorem I2.1[ b). In order to prove that the solution is 1-admissible, recall that 
c = " 2 ^" 2 ~^I^" 1 ~~ 1 ^ ■ Then, due to (14. ip . we need to show: 

v 2 (u 2 - 1) - Vi(Ui - 1) 
v 2 - v x 



Ai(«2, v 2 ) = u 2 - 1/2 - J 1/1 + v\ < 



< «1 - 1/2 - y/l/A + vl = Ai(tti,Vi). 

Since % = u 2 = w and i>i = 0, the latter reduces to 

u - 1/2 - yJl/4 + v% <u-l<u-l=^l/2- Jl/A + v~l < 0, 

which is clearly true. This concludes the proof. □ 

With this lemma established, we can attempt the proof of the following theorem. 

Theorem 4.1. Given any Riemann initial data (12.51) suc/i that v 2 < < v±, there 
exists a solution of (II. ip m t/ie sense of Definition \2. 1\ which consists of a combination 
of the classical Lax admissible simple waves (shock or rarefaction) and compressive 5 
waves, 1-admissible in the sense of Definition ^. 1\ 

Proof. The solution is plotted on Figure |2j First, we have the rarefaction wave-1 
(RW1) issuing from the left state L = (ui,Vi) and connecting it to the state (w m ,0). 
Then, we connect the state (u m , 0) with the state (u m ,v m ) by the 5-shock wave, 
and finally we connect (u m ,v m ) with R = (u 2 ,v 2 ) by the shock wave-2 (SW2) or 
rarefaction wave-2 (RW2). 

The solution is admissible since all the simple shocks which it contains are admis- 
sible. Namely, Lemma 14.11 provides admissibility for the 5-shock wave while other 
waves are admissible according to the standard theory (see Figured]). Furthermore, 
such combination of shocks is clearly possible since the speed of the state L equals 
to Xi(ui, v i) and it is less than the speed Xi(u m , 0) of the middle point (u m , 0) (since 
they are connected by the rarefaction wave). Furthermore, the speed of the 5 shock 
connecting (u m , 0) and (u m ,v m ) equals v m {u m — l)/v m = Xi(u m , 0) and it is slower 
than the speed of the state (u m ,v m ) which equals either X 2 ( y u m ,v m ) (if we have RW2 
between (u m ,v m ) and (u 2 ,v 2 )) or we have c = " 2 (" 2 ~ 1 jj^(" m ~ 1 ) < Ai(u m , 0) (if we 
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have SW2 between (u m ,v m ) and (u 2 ,v 2 )). Both situations are plotted in the phase 
space picture shown in Figure 

Indeed, if (u m , v m ) is connected to R = (112, v 2 ) by the RW2, then the speed of 
(u m ,v m ) is X 2 (u m ,v m ) > \\(u m ,0). On the other hand, if (u m ,v m ) is connected to 
R = (u2,V2) by the SW2, then its speed is 

v 2 (u 2 - 1) - v m (u m - 1) u 2 -u m 

= u m - 1 + t>2 > Ai(u m , 0) = u m - 1, 

^2 ~V m V 2 - V m 

since t>2 < 0, v 2 — v m > 0, and u 2 < u m (see Figure [2]). □ 




Figure 2. Thick full line: L ™ 1 (u m ,0) A (u m ,v m ) S ^ 2 R. Thick 
dashed line: L — >• (« m , 0) — > (u m , v m ) — > R. 

This theorem provides existence of an admissible 5-shock solution of the system 
( II. ip with Riemann data (12. 5p . However, even with the admissibility concept provided 
by Definition 14.11 it is not difficult to see that uniqueness may not hold. For example, 
a left state L = (u\,vi) and a right state R = (u 2 ,v 2 ) may be joined directly by a 
1-admissible 5 shock as long as 

u 2 - Ui 1 /T~ ~ u 2 - Ui 1 /T" ~ 

f 2 — f 1 2 V 4 v 2 —Vi 2 V 4 

and this is true whenever u\—u 2 is large enough and v 2 < < ui. We could, of course, 
add certain conditions which would eliminate the non-uniqueness. For instance, we 
could announce a S shock as admissible only if it connects states L = (u,0) and 
R = (u,v), v < 0. However, we do not have any physical justification for such 
condition and we shall confine ourselves on the existence statement. 
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Figure 3. The left plot corresponds to the situation when (u m ,v m ) is 
connected with {u2,v 2 ) by SW2, and the right plot when (u m ,v m ) is 
connected with {u 2 , v 2 ) by RW2. 



5. Appendix 

We shall end the paper by considering the possibility of the 5 distribution to be 
adjoint to the unknown v in (11.11) . We start with a lemma which will help us to 
connect certain states by admissible S shocks residing in the unknown v. 

Lemma 5.1. Assume that in (12.51) we have U\ = u 2 and V\ = —v 2 > 0. Then, if 
c = Xi(ui, v 1) = Aj(w 2 , v 2 ), i — 1, 2, the functions 

u(x, t) =Uq(x — ct), 

(5.1) 

v (x, t) =Vo(x — ct) + a(t)5(x — ct), 
where a(t) is given by (13. 5p . represent the i- admissible S-shock solution to (11.11) . 

Proof. It is enough to rely on Corollary [T] and proof of Theorem 12.11 Indeed, taking 
(15. ip for any c£l, and inserting them into Definition 14.1} we see that such u and 
v represent the 5-shock solution to (11.11) . (12.51) . To see this, one may use the same 
reasoning as in the proof of Theorem 12.11 and relation (I3.14p . 

Next, we take c = Ai(wi,t>i) = Ai(u 2 ,f2) or c = X 2 (ui,v 1 ) = \ 2 {u 2 ,v 2 ) to con- 
clude that the pair {u, v) is 1-admissible or 2-admissible, respectively, in the sense of 
Definition ED □ 

Using the lemma, we can connect the states 



L = {u\,vi) and R = {u 2 ,v 2 ) where u 2 > u\ 
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by an admissible 5-shock solution (u,v) to (II. ip . (I2.5P admitting the 5 shock in the 
function v through one of the following procedures 

• L — > (v m ,u m ) by RWI; (v m ,u m ) — > (—v m ,u m ) by the S shock with the speed 
c = Xi(u m ,v m ); (—v m ,u m ) — > R by RW2. See Figure HI In this case, we 
can also put c = \2(u m , v m ). If we take such c then 5 shock travels with the 
state (u m , —v m ). If we take c = Ai(w m , v m ) then 5 shock travels with the state 
(u m ,t> m ). Remark the non-uniqueness that we have here. 

• L — > (v m ,u m ) by SWI; (v m ,u m ) — > (—v m ,u m ) by the 5 shock with the speed 
c = X 2 (u m , v m ); (-v m ,u m ) -)■ R by RW2. See Figure EJ 

• L — )■ (v m ,u m ) by RWI; (v m ,u m ) — >• (— v m ,u m ) by the 5 shock with the speed 
c = Xi(u m ,v m ); (-v m ,u m ) -)■ i? by SW2. See Figured 




Figure 4. L = -> (w m , tf m ) ->• (w TO , -v m ) -> i?=(n 2 ,w 2 )- 

In the case when u 2 < Ui, we do not have a general recipe for connecting the 
states L = (ui,Vx) and R = (^2,^2) by an admissible 5 shock solution with the S 
function adjoined to v. Finally, observe that each of the 5 shocks in this section is 
not really compressive since characteristics from both sides of the shock are parallel 
to the shock. Thus we cannot say that concentration effects are present. 
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Figure 5. L = ->■ (w m ,v m ) ->■ (u m , -t> m ) ->■ it = {u 2 ,v 2 ) 
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Figure 6. L = (ui,vi) R ^r l (u m ,v m ) A (u m , -v m ) S ^ 2 R = (u 2 ,v 2 ). 
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